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SUMMARY 


A tension fatigue life prediction methodology for composite laminates is 
presented. Tension fatigue tests were conducted on quasi- i sotropi c and 
ortho tropic glass epoxy, graphite epoxy, and glass/graphite epoxy hybrid 
laminates. Edge delamination onset data were used to generate plots of strain 
energy release rate as a function of cycles to delamination onset. These plots 
were then used along with strain energy release rate analyses of delaminations 
initiating at matrix cracks to predict local delamination onset. Stiffness loss 
was measured experimentally to account for the accumulation of matrix cracks and 
for delamination growth. Fatigue failure was predicted by comparing the increase 
in global strain resulting from stiffness loss to the decrease in laminate 
failure strain resulting from delaminations forming at matrix cracks through the 
laminate thickness. Good agreement between measured and predicted lives 
indicated that the through-thi ckness damage accumulation model can accurately 
describe fatigue failure for laminates where the delamination onset behavior in 
fatigue is well characterized, and stiffness loss can be monitored in real time 
to account for damage growth. 


Keywords: Composite Materials, Matrix Cracking, Delamination, Fatigue, Stiffness 
Loss, Strain Energy Release Rate, Fracture Mechanics 
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Coefficient in power law for delamination growth 

Delamination size 

Laminate half- width 

uricracked ply thickness 

cracked ply thickness 

Axial modulus of a laminate 

Axial modulus before delamination 

Modulus of an edge delaminated laminate 
Modulus of a locally delaminated cross section 

Modulus of local cross section with edge and local delaminations 
initial modulus 

Lamina modulus in the fiber direction 

Lamina modulus transverse to the fiber direction 

In-plane shear modulus 

Strain energy release rate 
Mode I strain energy release rate 

Mode II strain energy release rate 

Critical value of G at delamination onset 

Maximum G in fatigue cycle 

Strain concentration factor 

Laminate length 

Number of sublaminates formed by edge delamination 
olope of G versus log N curve for delamination onset 
Exponent in power law for delamination growth 
Number of fatigue cycles 
Cycles at failure in fatigue 

Number of local delaminations through the laminate thickness 
Number of load steps 

Cyclic stress ratio in fatigue (a . /a ) 

min max 

Matrix crack spacing 
Thickness 

Laminate thickness 

Thickness of a locally delaminated cross section 
Uniaxial strain 

Critical strain at delamination onset 
Strain at failure 
Maximum strain in fatigue cycle 
Uniaxial stress 

Maximum stress in fatigue cycle 
Minimum stress in fatigue cycle 

2 


1 



ORIGINAL PAGE IS 
OF POOR QUAUTY 


INTRODUCTION 

As composite materials are considered for primary structural applications, 
the need to predict their fatigue behavior becomes important* For continuous 
fiber reinforced composites, delamination is the most commonly observed damage 
mechanism. Fracture mechanics analyses of delaminations have been shown to 
provide generic characteri zat ions of delamination onset and growth [1-4]. 
However, fatigue life prediction requires not only a character i zat ion of 
delamination behavior, but an assessment of the consequence of fatigue damage as 
well [ 3 “8 3 - The purpose of this paper fs to develop a fatigue life prediction 
methodology for composite laminates subjected to tension fatigue loading that 
incorporates both the generic fracture mechanics char act eri zat ion of 
delamination as well as the assessment of the influence of damage on laminate 
fatigue life. 

MATERIALS 

Composite laminates were made from pre-impregnated rolls of T190 E~glass fibers 

embedded in an X751/50 epoxy matrix, and HM35 graphite fibers embedded in the 

same X751/50 epoxy matrix. The prepreg material was supplied by Cyanamid 

Fothergill, Ltd., U.K.. Laminates were cured in a press by Agusta Helicopter 

Co., Italy, using the manufacturer’s recommended curing cycle consisting of a 

2.7 atm pressure and a maximum temperature of 125°C. Glass epoxy, graphite 

epoxy, and glass/graphite epoxy hybrid laminates were manufactured with an 

average glass ply thickness of 0.226 mm, graphite ply thickness of 0.1875 mm, 

and a 48-53$ fiber volume fraction. Specimens 38.1 mm wide by 300 mm long were 

cut from the plates after they were cured and were subsequently tested with a 

180 mm length between the grips of the load frame. Tests were conducted on [0] Q , 

o 

[90]^2> anc * [45/~45 ]^ s glass epoxy and graphite epoxy laminates with similar ply 

thicknesses and fiber volume fractions to generate the initial lamina property 
data shown in tables 1-3* Two different batches of T190 E-Glass X751 /50 epoxy 
prepreg were supplied. The first batch was used to make the glass epoxy 
laminates (table 1), and the second batch was used to make the hybrid laminates 
(table 2). Although the lamina moduli were similar for these two batches, the 
zero degree failure strains were significantly lower for the first batch (tables 
1 & 2 ) . 

EXPERIMENTS 

Static and fatigue tests were performed on ( 45/-45/0/90) s quasi-isotropic, 

glass epoxy and graphite epoxy laminates. In addition, static and fatigue tests 
were performed on two hybrid laminates. The first was a ( 45/-45/0/90) quasi- 

isotropic laminate where the 45, ”45, and 0 degree plies had graphite fibers and 
the 90 degree plies had glass fibers in the same X751/50 epoxy matrix. The 
second was a (4b/-45/0) orthotropic laminate where the 45 and -45 degree plies 

were graphite and the 0 degree plies were glass with the same X751 /50 epoxy 
matrix . 

Static tension and t ens i on- t ens i on fatigue tests were conducted on an 
Instron ser vo-hydr aul i c test stand. For the static tests, load was applied 
axially in stroke control. Fatigue tests were conducted in load control by 
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applying a sinusoidal load about the mean load at a frequency of 6 Hz and an r — 
ratio of 0.1. All tests were conducted in an ambient environment with the 
temperature ranging between 18° and 25°C. Axial strain was derived from 
displacements measured by a pair of Direct Current Differential Transducers 
(DCDT's) over an 82 mm gage length. Output signals from the load cell and the 
DCDT’s were plotted on an X-Y recorder. The initial slope of these load 
d i spl acemen t plots were used to calculate laminate moduli. During the static 
tests the machine was stopped at fixed intervals and the damage was recorded. 
For the glass epoxy laminates, a high intensity light was shone through the 
specimen to observe and photograph the damage. Because the X751 epoxy was 
opaque, a colored dye penetrant was applied to the specimen edges to provide a 
color-enhanced contrast for viewing the damage. Similarly, a zinc-iodide dye 
penetrant was used to enhance the X-ray images taken of the graphite epoxy and 
hybrid laminates for viewing the damage. Matrix crack density and delamination 
size were measured from high contrast photographic prints. The delaminated area 
was measured using a graphic tablet to trace the shape of the delamination. This 
information was transferred to the computer software to calculate the enclosed 
area . 

For each laminate, five quasi-static tension tests were run with a ram 
speed of 0.5 mm/min. At specified loads, the specimens were unloaded to record 
damage and residual stiffness. During the tension— tension fatigue tests, cyclic 
loading was interrupted at specified intervals to monitor damage and measure 
static stiffness. To evaluate stiffness, load versus displacement curves were 
generated by applying a sinusoidal cycle at 0.006 Hz. 

STATIC TENSION BEHAVIOR 

For the quasi-isotropic graphite epoxy laminates, edge delaminations formed 
in the 0/90 interfaces before final failure of the laminate. Matrix cracks were 
also visible in the 90 degree plies, however, they had no noticeable effect on 
the laminate modulus before edge delamination onset. The stress-strai n curves 
were linear up to the onset of edge delamination (fig.1). 

For the quasi-isotropic glass epoxy laminates, edge delaminations also 
formed in the 0/90 interfaces before final failure of the laminate. Matrix 
cracks were also visible in the 90 degree glass plies. However, unlike the 
graphite epoxy laminates, the stress-strain curves became nonlinear after the 
onset of 90 degree matrix cracking (fig.1). 

For the qua a i - i sot r op i c hybrid laminates, static tests showed no 
edge delamination adore final fracture of the laminate. Matrix cracks were 
visible in the 90 degree glass plies, however, they had no noticeable effect on 
the laminate modulus. The stress-strai n curves were linear up to the final 
failure (fig. 2). The final failure strains for the quas i - i sot ropi c hybrid 
laminates (table 4> were similar to the failure strains for the high modulus 
graphite fiDers (table 3). 

For the orthotropic hybrid laminates, no matrix cracking or delamination 
was observed in the static tests before final failure of the laminates. However, 
the stress-strain curves were nonlinear (fig. 2). The average failure strain for 
the orthotropic laminates (table *4) was similar to the failure strain for the E- 
glass fibers (table 2). 

TENSION FATIGUE BEHAVIOR 

Figure 3 shows three dye penetrant enhanced photographs of the tension 
fatigue damage in the ( 45/-45/0/90) glass epoxy laminates taken at different 
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stages of cyclic loading. These photographs clearly show the edge delaminations 
and matrix cracks in the 45 degree plies. Local delaminations can also be seen 
extending from the 45 degree matrix cracks. Figure 4 shows a schematic of some 
of this damage, including edge delaminations that form at the edge in the 0/90 
interface and jump through 90 degree ply cracks to the other 0/90 interface, and 
local delaminations that form in the 45/-45 interface, originating at 45 degree 
matrix ply cracks . 

The fatigue damage in the glass epoxy laminates progressed in the following 
sequence. First, extensive matrix cracking developed in the 90 degree plies, 
followed by edge delamination in the 0/90 interfaces. Next, matrix cracks 
appeared in the 45 degree and -45 degree plies and initiated local 
delaminations, first in the 45/~45 interfaces, followed by the -45/0 interfaces. 
Finally, after enough local delaminations had formed through the thickness at a 
particular location, fiber failure occurred and the laminate failed. The same 
damage mechanisms, and similar sequences of damage, were observed in the 
graphite epoxy and hybrid quasi-isotropic laminates. Dye penetrant enhanced 
radiographs of this damage in T300/5208 graphite epoxy laminates have been 
documented previously [6,7]. 

For the glass epoxy laminates, fig. 5 shows the number of cycles at a given 
maximum cyclic stress for edge delaminations to form (solid symbols), for the 
first local delamination to form at the 45/-45 interface (brackets), and for 
fatigue failure to occur (open symbols). Fatigue failures occurred in les3 than 

10^ cycles at maximum stress levels on the order of 50? of the static strength 
or greater. 

For the graphite epoxy quasi-isotropic laminates, fig. 6 shows the static 
ultimate strength and the maximum stress versus cycles to onset of edge 

delamination. Edge de 1 am i nat i ons formed in less than 10^ cycles at maximum 
stress levels above 50? of the static strength, but no fatigue failures occurred 

before 10^ cycles. Three fatigue tests run at 90? of the static ultimate stress, 
where edge delaminations occurred in the first load cycle, also had fatigue 

0 

lives greater than 10 cycles. Although the entire sequence of failure in 
graphite epoxy laminates has been studied in the past [6,7], for this study only 
the data up to delamination onset was used for characterization of delamination 
onset in fatigue. 

For the quasi-isotropic, hybrid laminates, fig. 7 shows the maximum cyclic 
stress as a function of cycles to delami nation onset and cycles to failure. 
Although edge delaminations did not form in the hybrid quasi-isotropic laminates 
under static loading, they did form during cyclic loading at stress levels below 

the static strengtn of the laminate. Edge delaminations formed in less than 10^ 
cycles at maximum stress levels above 50$ of the static strength. However, 

fatigue failures at less than 10^ cycles were only observed at maximum stress 
levels greater than 90$ of the static strength. At these high stress levels, 
local delaminations also occurred in the 45/~45 interface, growing from matrix 
cracks in the outermost 45 degree plies. Although these local delaminations were 
visible on the specimen during testing, they were difficult to document using 
the dye penetrant enhanced radiography because the small local delaminations 
were obscured in the radiographs by the image of the 0/90 interface edge 
delaminations . 

For the (45/-45/0) s hybrid laminates, fig. 8 shows the maximum cyclic stress 

as a function of cycles to onset of (1) matrix cracks in the 45 or -45 degree 
plies, (2) the onset of local delamintion in the 45/-45 degree interface, and 
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(3) the final fatigue failure. Fatigue failures occurred in less than 10^ cycles 
at maximum stress levels on the order of 50 % of the static strength or greater. 
These fatigue failures were precipitated by the onset of local delaminations in 
the 45/-45 interface that grew from matrix cracks in the 45 and -45 degree 
plies. Sections of the failed composite were polished and examined in an optical 
microscope to verify the presence of these damage mechanisms. These sections 
also showed that local delaminations had formed in the -45/0 interfaces. These 
local delaminations were not visible on the radiographs. Furthermore, these 
polished sections were used to document whether the delamination and matrix 
cracking was occurring in the matrix or at the fiber matrix interfaces. In the 
graphite plies, matrix cracks occurred primarily within the matrix, and only 
occasionally at the fiber matrix interface. However, in the glass epoxy 
laminates, and in the hybrid laminates, the fractographs showed that the 
delaminations occurred almost exclusively within the glass plies at the fiber 
matrix interfaces closest to the ply interface (fig. 9). 

In order to predict the fatigue failure of these laminates, the onset and 
growth of the damage observed must be characterized , and the influence of this 
damage on laminate stiffness and strength must be determined. To this end, a 
fatigue analysis and life prediction methodology. 

FATIGUE ANALYSIS AND LIFE PREDICTION METHODOLOGY 


Influence of Damage on Laminate Stiffness 


Several damage mechanisms may reduce the stiffness of a composite laminate. 
The most commonly observed fatigue damage mechanisms are matrix cracking and 
delamination. In this section, the influence of matrix cracking and delamination 
on laminate stiffness will be reviewed. 

Fig. 10 shows the influence of damage on laminate stiffness. As matrix 
cracks accumul ate , and as delaminations form and grow, the stiffness of the 
laminate decreases. Laminate stiffness is the ratio of the remote stress to the 
global strain in the laminate. This global strain is typically measured using an 
extensometer or DC!)T (fig. 10), which yields the displacement of the laminate 
over a fairly long gage length relative to the laminate's length. As damage 
forms and grows in the laminate under a constant maximum cyclic stress, 
correspondi ng to a constant applied maximum cyclic load, the global strain in 
the laminate increases. 


Previous studies have determined the relationships between stiffness loss 
and damage extent [1-10J. The amount of stiffness loss associated with matrix 
cracking depends upon the ply orientation of the cracked ply, the laminate 
layup, the relative moduli of the fiber and the matrix, and the crack spacing, 
or density of cracks, in the ply. For example, in ref. 4, an equation was derived 
for stiffness loss due to matrix cracking in 90 degree plies of cross ply 
laminates as 


e lam 
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and c and d are the thickness of the cracked and uncracked plies, respect ively . 
As the crack density increases, i.e., as the crack spacing, 2s, decreases, the 
stiffness of the laminate will decrease. For a particular layup with a given 
crack density, the? stiffness los3 would be larger for’ glass epoxy than graphite 
epoxy because the ratio of the transverse to axial modulus, ^22^1 1 ’ is greater 

for glass epoxy than for graphite epoxy (fig. 11a). For example, fig. 11b shows 
the maximum stiffness loss when matrix cracking has saturated in the 90 degree 
plies of cross ply laminates made of graphite epoxy and E-glass epoxy. The glass 
epoxy laminates have significantly greater stiffness loss than the graphite 
epoxy laminates. 

The amount of stiffness loss due to delamination also depends on the 
laminate layup and the relative modulus of the fiber and the matrix, as well as 
the location and extent of the delamination. As delaminations form and grow in a 
particular interface, the laminate stiffness decreases as the delamination size, 
a, increases. In ref.1, an equation was derived for the stiffness loss 
associated with edge delamination as 

E = ( E*~ E t )a/b + E. (3) 
LAM LAM 


where a/b is the ratio of the delamination size to the laminate half-width, and 

* 

E is determined from a rule of mixtures expression 


# 
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( 4 ) 


where the moduli, E. , of the M sublaminates formed by the delamination are 

l 

# 

calculated from laminated plate theory. The difference in E and E reflects 

the loss of transverse constraint in the sublaminates formed by the 
delamination. The greater the ratio of the transverse to axial lamina moduli, 
E 22 /E 11* ies3 transverse constraint, resulting in less difference in 

* 

E. and E . For example, fig. 12 shows the percentage stiffness loss, as 
LAM 

# 

measured by (E f - E )/E t „ ^ , for edge delamination in the 0/90 interfaces of 
LAM LAM 

(45/-45/0/90) E-glass epoxy, graphite epoxy, and hybrid laminates. As noted 

earlier, the glass epoxy has the largest ratio of anci hence ^ as 

smallest stiffness loss when delaminated. 

Delaminations starting from matrix cracks will effect laminate stiffness 
differently than delaminations growing from the straight edge. In ref. 2, an 
equation was derived for the stiffness loss associated with delaminations from 
matrix cracks as 


E 


{ (a/£)t LANJ Cl / ( t LE) E LD ) 


1 / ( t. AM E t ) 
LAM LAM 


] } 


-1 


( 5 ) 


where a/l is the ratio of the delamination length to the laminate length, and 

t and E represent the thickness and modulus of the locally delaminated 
L D L D 
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region in the vicinity of the matrix crack. Hence, the locally delaminated 
modulus, E ld , is the modulus of the laminate without the cracked plies that 

caused the local delamination. This locally delaminated modulus is similar to E* 
in eq(4). However, in addition to reflecting the loss in transverse constraint 
due to the delamination, E LD also reflects the loss of the load bearing capacity 

of the cracked ply. Similar to edge delamination, the st ffness of the laminate 
decreases as the size of the delamination increases. However, unlike edge 
delaminations which form at the two edges and grow progressi vely towards the 
center of the laminate width, local delaminations tend to accumulate at several 
matrix cracks along the length, growing only a smal' distance at any one 
location. The cumulative effect of these local delam nations with cycles, 
however, may have a significant effect on measured stiffness loss. 

Delamination Onset and Growth Characteri zation 


In order to predict stiffness loss as a function of fatigue cycles, the 
onset and growth of matrix cracks and delaminations must be characterized in 
terms of a generic parameter that is representative of the composite material 
being tested, but independent of laminate structural variables such as layup, 
stacking sequence, and ply thickness. This characterization is typically 
achieved using the strain energy release rate, G, associated with matrix 
cracking and delamination [1-4,11-17]. 

Fig. 13 shows the steps that would be required to predict stiffness loss as 
a function of fatigue cycles using a G character i zation of damage onset and 
growth. First, plots of G versus log N must be generated to characterize the 
onset of delamination [3,12-16], and power law relationships between G and the 
rate of growth of delamination with fatigue cycles are needed to characterize 
damage growth [11,15] (fig.1 33 ). Similar appr oaches may be used t o char act er i ze 
the onset and growth of matrix cracks [4]. Using these material 
characterizations, the decrease in matrix crack spacing, 2s, and the increase in 
delamination size, a, with fatigue cycles may be predicted (fig. 13b). This 
information, in turn, may be used with eqs(1-5) to predict the decrease in 
modulus with cycles, which for a constant stress amplitude test is tantamount to 
predicting the increase in global strain with cycles (fig. 13 c). 

Although this technique may be demonstrated for cases where there is one 
dominant damage mechanism, application of this approach in general is difficult 
because the various damage modes interact, complicating their unique 
characterization in terms of G [4]. For example, although the elastic analysis 
for G associated with edge delamination growth is independent of delamination 
size, stable delamination growth is often observed experimentally [1,6,17]. The 
strain energy release rate for edge delamination is given by 


G 


l S- tH ( E - E *) 

2 ' LAM j 


( 6 ) 


which is independent of the delamination size. Theoretically, when a critical 

value of strain, e , is reached, c or r es po nd i ng to a critical G . the 

c c 
delamination should form on the edge and grow immediately through the width. 
However, edge de 1 am i nat i ons usually grow in a stable fashion, requiring 
increasing 3train levels, and hence increasing G, for the delamination to grow 
across the width. This stable growth may be correlated with the accumulation of 
90 degree matrix cracks ahead of the delamination front. For example, fig. 14 
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shows a plot of normalized delamination size, a/b, as a function of the strain 
applied to an eleven-ply ( ±30/±30/90/90) T300/5208 graphite epoxy laminate. The 

3 

edge delaminations form at a strain of approximately 0.0035, but do not grow 
across the specimen width until the strain reaches approximately 0.0065. Also 
plotted in fig. 1*1 on the right hand ordinate is the 90 degree ply crack spacing 
measured in the center of the laminate. There appears to be a direct correlation 
between the stable delamination growth and the accumulation of matrix cracks 
ahead of the delamination front. These matrix cracks apparently interact with 
the delamination and increase the resistance to delamination growth. Stable 
delamination growth may be predicted by generating a delamination resistance Fl- 
eur ve using eq.6 [1,17]. However, the resulting R-curve is no longer generic, 
because the matrix cracking that is causing the delamination resistance is 
governed by structural variables such as ply thickness and stacking sequence. 

Delamination also influences the formation and accumulation of matrix 
cracks. Delamination relaxes the constraint of neighboring plies, and hence 
changes the saturation spacing of matrix cracks in the of f-axi s_ pi ies . For 
example, when delaminations form at the edges of the (±30/±30/90/90) s laminate 

shown in fig. 6, the constraint between the -30 deg and 90 deg plies is relaxed, 
and the 90 deg cracks form sooner, with smaller crack spacings than possible if 
no delamination had existed [1]. An R-curve description of matrix cracking has 
been used to describe the accumulation of matrix cracks, similar to the approach 
that has been attempted for delamination [*J]. However, when these cracks 
interact with delaminations, this description is no longer generic. 

Even if one could achieve a truly generic description of damage 
accumulation with cycles, the resulting stiffness loss prediction, and hence the 
prediction of increasing global strain with cycles, is necessary, but not 
sufficient, to predict fatigue life. The final failure of the laminate is 
governed not only by loss in stiffness, but also by the the local strain 
concentrations that develop in the primary load bearing plies, which in most 
laminates are zero degree plies. 

Influence of Local Strain Concentrations on Failure 


Fig. 15a shows that fatigue failures typically occur after the global strain 

has increased because of the fatigue damage growth, but before this global 

strain reaches the global strain at failure, e , measured during a static 

r 

strength test [5-3]. Therefore, local strain concentrations must be present in 
the zero degree plies that control the laminate strength. Although matrix cracks 
create small strain concentrations in the neighboring plies, their magnitudes 
are generally small because the stiffness of the cracked ply is usually much 
less than the stiffness of the zero degree ply [8]. Furthermore, strain 
concentrations due to matrix cracks act over only a local volume near the crack 
tip in the adjacent pLy [8J. Hence, the final failure in a zero degree ply of a 
laminate may follow a neighboring ply crack [18], but the laminate failure 
strain will not be strongly influenced by the presence of the matrix cracking 
[8J. Once delaminations initiate at matrix ply cracks, however, the local strain 
will increase significantly throughout the remaining through-thickness cross 
section [2,5-8]. These local strain increases may not be sensed by the global 
strain measurement, because delaminations starting from matrix cracks grow very 
little once they form. If several delaminations form at matrix cracks throughout 
the laminate thickness at one location, then the local strain on the zero degree 
plies at that location may reach the static failure strain, resulting in the 
observed fatigue failure (fig. 15b). 
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This mechanism for fatigue failure has been observed previously for 
graphite epoxy laminates [3,5-8], and the local strain concentrations resulting 
from cumulative local delaminations through the thickness have been quantified 
[2, 3,5-8]. These local strain concentrations may be calculated as 


K 

e 


E LAM t LAM 

E LD t LD 


( 7 ) 


Typically, the local strain concentration will result in a trade off between the 
increased modulus, E^ D > E LAM> which occurs because E^ D is a more zero degree 

dominated layup than the original laminate, and the decrease in load bearing 
cross section, t^< t . Fig. 16 shows the local strain concentrations 

calculated for local delaminations accumulating through the thickness in a 
( ^45/— ^5/0 ) s E-glass epoxy laminate. For each progressive delamination that forms 

through the thickness, the local strain concentration on the zero degree plies 
increases. 

In some cases, these local delaminations will be influenced by other 
mechanisms. For example, if the local delamination forms at an edge, where an 
edge delamination exists in another interface (fig. 17), the local strain 
concentration in the zero degree plies near the edge will be 


„* E t LAM 

l\ W 

^LD^LD 


( 8 ) 


instead of calculated from eq.7, which is valid away from the edge 
delamination in the interior of the laminate width. Fig. 18 shows values of K 

e 

* 

and for thr ough-the-thi ckness accumulation of local delaminations in 

( 4b/-45/0/90 ) E-glass epoxy laminates containing delaminations near the edges 

in the 0/90 interfaces. For a given through-thi ckness accumulation of local 

delaminations, < K _ , because E /E^ D < E LAM /E LD - 1° fig-19, K* values are 

plotted for the same quas i - i sot r opi c laminate using both E-glass epoxy and 
graphite epoxy properties. As shown in fig. 19, for the same layup and damage 
situation, the resulting local strain concentration will be more severe for the 
E-glass epoxy than for the graphite epoxy laminate. 

As noted earlier, each time a local delamination initiates from a matrix 
crack, the local strain in the the remaining through thickness cross section, 

and hence in the zero degree plies, increases by an amount equal to K times the 

e 

global cyclic strain, e , until the local strain reaches the static failure 

m a X 

strain, (fig. 20a). A simpler way to visualize this process however, is to 
reduce the static failure strain to some effective global e value each time a 

r 

new local delamination forms through the thickness. Hence, the effective failure 

strain for the "i"th delamination, (e )., would be equal to e„/(K ).. As local 

r i F e i 

delaminations accumulated through the thickness, (e„). would decrease 

r 1 

incrementally. Fatigue failure would correspond to the number of cycles where 
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the damage growth increased the global maximum cyclic strain to the current 

value of (z ). (fig. 20b). This approach does not require a prediction of damage 
F i 

growth with fatigue cycles if the laminate stiffness loss, and hence the 

increase in global strain, can be monitored in real time. When this is possible, 

only the incremental decreases in (e_). needs to be predicted to predict fatigue 

F 1 

life. This may be accomplished by first assuming that matrix cracks exist in all 
of the off-axis plies. Then, the number of fatigue cycles to onset of each local 
delamination through the thickness may be predicted using delamination onset 
criteria (fig. 13 a) along with strain energy release rate analyses for local 
delamination. At a given location, as each local delamination forms through the 
thickness, may be reduced by the appropriate and compared to the current 

value of z , based on measured stiffness loss, to determine if fatigue failure 
max 

will occur. Hence, the ability to predict local delamination onset, and its 

effect on facilitates using measured stiffness loss to predict fatigue life. 
F 

Tension Fatigue Life Prediction Procedure 


In order to predict the onset of local delaminations with fatigue cycles, 
the G versus log N characterization of the composite material must be 
generated. Data from several materials with brittle and tough matrices indicate 

r< £ 

that between 10 < N < 10 cycles, the maximum cyclic G may be represented as a 
linear function of log N (fig.21), where N is the number of cycles to 
delamination onset at a prescribed G max L 1 3 1 - Hence, 

G = m log N + G (9) 

c 


where G and m are material parameters that characterize the onset of 

delamination under static and cyclic loading in the material (fig. 13a). This 
characterization may be accomplished using a variety of interlaminar fracture 
test methods L12-16J. Next, G must be calculated for each local delamination as 
it accumulates through the thickness. The first local delamination will 
typically form at a matrix crack in the surface ply, but the appropriate 
sequence may be confirmed by calculating G for matrix cracking in all of the 
off-axis plies in the laminate. The one with the highest G for the same applied 
load will be the first to form. This G may be calculated using the equation 
derived in ref. 2 for the strain energy release rate associated with local 
delaminations initiating at matrix cracks 


2 fc 2 

i/U ld e ld ) - ,/( Wlam ) j 


no) 


As shown in fig. 22, the thickness and modulus terms in eq(10) change for each 
succesive local delamination that forms through the thickness. For example, t LD 

and E for a 45/-45 local delamination in a (45/-45/0) laminate becomes the 
L D s 

t and E values used to calculate G for the next local delamination that 
LAM LAM 

forms through the thickness. Therefore, under a constant o , as local 

IllaX 

delaminations accumulate through the thickness, the driving force (i.e., G) for 
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each new delamination changes. Hence, fatigue life prediction for composite 
laminates requires a "cumulative damage" calculation, even for constant 
amplitude loading. To calculate the number of cycles for each succesive local 
delamination to form, N. , the appropriate expression for G from eq(10) is set 

equal to the delamination onset criterion of eq(9) (fig. 22) and the equation is 
solved for (fig. 23). Hence, 


log N, . i I V* l/(t LD E LD ) - ]. - G o l (,,) 

and the effective (e_). = e_/(K ).. 

F 1 F e l 

Fatigue failure will occur when e m . )X - (e^h, resulting in a fatigue life, N , 
of 


n f - 


N. 


( 12 ) 


i = 1 


where p is the number of local delaminations that form through the thickness of 
the laminate before failure. 

Because of the scatter in experimental data, the constant load amplitude 
fatigue life prediction methodology outlined in fig. 23 more closely resembles 

fig. 2^1. The variation in initial laminate modulus (i.e., the variation in c ) 

max 

and the variations in static failure strains from specimen to specimen must be 

taken into account. Hence, a range of possible fatigue lives would be predicted, 

rather than a single value. The lowest life would occur when then minimum value 

of (e ) in the e distribution reaches the largest e value in the 

p 1 f max 

distribution resulting from variations in laminate moduli. 


KE5ULTS 


Fatigue Delamination Characterization 


The fatigue life methodology outlined in the previous sections was used to 
predict the fatigue life of the quasi-isotropic glass epoxy laminates and the 
quasi-isotropic and orthotropic hybrid laminates. To make these predictions, the 
delamination onset behavior in fatigue for these materials was first charac- 
terized in terms of strain energy release rates. 

For the glass epoxy quasi-isotropic laminates, the lamina properties used 

* 

to calculate E LftM and E were reduced to reflect the stiffness loss that 

preceeded edge delamination onset. Specifically, E^ 2 was reduced by 29% based on 

a shear lag analysis [1,9] to reflect the saturation cracking in the 90 deg 
plies. However, this reduction in E 22 of the 90 degree plies, when substituted 

into the laminated plate theory, yielded laminate moduli that were significantly 
greater than moduli measured just before edge delamination onset. Because the 
±4b degree plies were loaded to relatively high strain levels (>0.01) before 
edge delamination onset, the nonlinearity in the stress-strain curve before edge 
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delamination onset was attributed to the matrix dominated response of the ±45 
degree glass epoxy plies. This nonlinearity is similar to the nonlinearity 
observed in quasi-static tests of (+45/0) E-glass X751/50 epoxy laminates 

3 

which show no evidence of matrix cracking in the ±45 degree plies. Hence, in 

the ±45 degree pLies was reduced by b>0% to reflect the nonlinearity observed. 
These reduced lamina properties (table 1) yielded quasi-isotropic longitudinal 
moduli, calculated by laminated plate theory, that were consistant with E^ 

values measured just before edge delamination onset in the fatigue tests (table 
4). In the absence of real time stiffness measurement, the cyclic plasticity of 
the ±45 degree plies would have to be char ac t er i ze d in addition to the 
characterization of matrix cracking and delamination [19]. However, this would 
only be necessary in laminates where the zero degree plies were loaded to very 
high cyclic strain levels without failing, as was the case for the :ero degree 
glass plies in the quasi-isotropic glass epoxy and orthotropic hybrid laminates 
( table 4) . 

The maximum cyclic strains for edge delamination onset in the quasi- 
isotropic glass epoxy laminates were used in eq(6) to calculate strain energy 
release rates, and the data were plotted versus the number of cycles to 
delamination onset (fig. 25). There was significant scatter in the static data 
for G c » possibly due to the interaction that occurred between the edge 

delamination as it formed, and the 90 deg ply cracks that were extensive before 
edge delamination onset. However, values from edge delamination data may be 

artificially elevated if extensive 90 deg ply cracking is present in the 
laminate [4,12]. Therefore, the minimum values in fatigue were used in eq(9) to 
characterize delamination onset and extrapolate an effective static G c . All of 

these data were generated at the 6 Hz frequency except for the N=2 data point. 
This data point was determined from a slow sinusoidal loading where the edge 
delamination formed at the second peak in the wave form. For the E-glass X751/50 

2 

epoxy, an effective G ^ value of 9 6.2 J/m was obtained, and the slope, m, was - 
10.52. 

Because the stress-strain behavior of the graphite epoxy quasi-isotropic 
laminates was linear up to the onset of edge delamination, the delaminated 

# 

modulus, E, could be obtained directly from the experimental data. Although 
90 degree ply cracks form after the onset of edge delamination and will 
contribute to stiffness loss, the influence of these matrix cracks is not large 
for graphite epoxy laminates (fig. 11b). However, this influence will be 

# 

reflected in the experimentally determined value for E . Fig. 26 shows the loss 

in stiffness as a function of delamination area for the gr aphi te/ epoxy 

laminates. Data from both the static and fatigue tests were plotted, and were 

correlated using a linear least squares regression fit. The extrapolation of the 

linear fit to total delamination, a/b*1.0 f yielded the delaminated modulus. 

Fig. 27 shows the maximum strain energy release rate as a function of cycles 

to edge delamination onset for the quasi-isotropic graphite epoxy laminates. As 

has been prevously noted [6,7,12], 90 degree ply cracking is not extensive in 

( 45/-45/0/90) graphite epoxy laminates before edge delamination onset. Hence, 
s 

the static G values for grahite epoxy have less scatter than was observed for 
c 

the glass epoxy laminates, and these G^ values are not artificially elevated. 
Hence, a linear relationship for eq,9 may be easily fit to the data using the 
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minimum static and fatigue delamination onset values measured (fig. 27) . For the 
HM35 graphite X751 /50 epoxy, a G q value of 137 J/m 2 was obtained, and the slope, 
m, was -10.5. 

Fig. 28 compares the maximum strain energy release rate as a function of 
cycles to edge delamination onset for the quasi-isotropic graphite epoxy and 
glass epoxy laminates. The slope, m, was the same for both materials. However, 
although both materials consisted of the same X751 epoxy matrix, the glass/epoxy 
curves are slightly lower than the graphite/epoxy ones. These lower values were 
attributed to the glass epoxy delaminating at the fiber-matrix interfaces within 
the 90 degree plies, instead of failing within the matrix layer between the zero 
degree and 90 degree plies, similar to the interfacial failure of the hybrid 
laminates as observed using fractogr^phy (fig. 9). 

Fig. 29 shows the maximum strain energy release rate as a function of cycles 
to edge delamination onset for the quasi-isotropic hybrid laminates. As noted 
earlier, edge delamination did not occur under static loading in these hybrid 
laminates. Hence, a linear fit to eq.9 was not possible. However, fig. 30 shows 
that the maximum strain energy release rate as a function of cycles to edge 
delamination onset was similar for the quasi-isotropic hybrid and glass epoxy 

laminates between 10 and 10 cycles. This was not surprising because the edge 
delaminations that formed in the 90 degree glass epoxy plies in both laminates 
were interfacial failures, as evidenced by fiber-matrix separation within the 90 

degree ply adjacent to the 0/90 ply interface. Hence, the effective G and m for 

c 

the hybrid quas i - i sotropic laminates was assumed to be identical to the values 
obtained for the glass epoxy laminates. 

Tension Fatigue Life Prediction 


To predict fatigue life for the quasi-isotropic glass epoxy laminates and 
the quasi-isotropic and orthotropic hybrid laminates, the onset of local 
delaminations from matrix cracks through the laminate thickness was first 
predicted using eq.11. Local delaminations that occurred between graphite epoxy 
plies were predicted using G c and m values for graphite epoxy; whereas, local 

delaminations that occurred between glass epoxy plies, or between adjacent 
graphite and glass epoxy plies, were predicted using the effective G and m 

G 

values for glass epoxy because they exhibited similar interfacial failures. For 
the quas i - i sot r op i c laminates, local delaminations formed at the edge and 
quickly grew normal to the matrix crack along the laminate width (fig. 17). 
Hence, for the quasi-isotropic laminates, E^ and E LD in eq.11 were replaced by 

# # 

E and E^, respectively, because local delaminations initiated at the edges 

after edge delaminations had already formed in the 0/90 interfaces. 

For the glass epoxy quasi-isotropic laminates, reduced values of E 22 and 

G 12 discussed previously (table 1) were used to calculate the moduli in eq(11). 


For the quas i - i sot ropi c hybrid laminates, original lamina properties (tables 
2&3) were used to calculate moduli in eq.11 because no stiffness loss was 
observed before delamination onset in these laminates. For the orthotropic 
laminates, the ±45 degree graphite epoxy plies developed matrix ply cracks that 
reduced the laminate modulus before they led to the formation of local 
delaminations. Furthermore, the non-linearity in the stress— strai n curves 
observed under static tests may have also affected the modulus measured in 
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fatigue. Hence, the lamina shear modulus in these layers was reduced by 50% 
(table 3) so that the (45/-45/0) s ^AM P re< ^^ c ^ e( ^ from laminated plate theory 

would agree with the modulus, E Q , measured in fatigue just before the onset of 

local delamination (table 4). 

Reductions in effective e for local delaminations accumulating through the 

r 

thickness were calculated using eq.7, wher^ E. and K fA were calculated from 

LAM Li) 

laminated plate theory. As previously noted, edge delaminations formed before 
local delaminations in the quasi-isotropic laminates. Therefore, initially the 
strain concentration for a local delamination starting at the edge of a quasi- 

# 

isotropic laminate is K , as determined by eq.8. However, this strain 

concentration quickly becomes when the local delamination forms in the 

interior. Hence, K values were used to calculate reductions in e r , for all the 

e P 

laminates. 

Figs. 31 "33 show the range of estimated and measured fatigue lives for the 
quasi-isotropic glass epoxy laminates, the quasi-isotropic hybrid laminates, and 
orthotropic hybrid laminates. In some cases the estimated lives are 
conservative, but in general, very good agreement was observed. 

DISCUSSION 

The good agreement between measured and predicted lives (figs. 31 “33) 
indicates that the through- thi ckness damage accumulation model can accurately 
describe fatigue failure for a material whose delamination behavior in fatigue 
is well char act er i zed , and if stiffness loss is monitored in real time to 
account for damage growth. 

For the glass epoxy laminates, the G vs. log char act er i zat i on was 
generated using data from the same laminates whose fatigue lives were being 
predicted. The fatigue lives of the hybrid laminates, however, were predicted 
using the G vs. log N oharacteri zation from other glass epoxy and graphite epoxy 
laminates having the same constituents (fibers and matrix) as the hybrids. 
Ideally, the G vs. log N characterization should be performed using 
standardized laboratory tests [15], and should then be used to predict the 
fatigue behavior of structural components made of the same material. Hence, for 
the fatigue life prediction to be accurate, the laboratory characterization 
should be performed on identical materials (same constituents, fiber volume 
fraction, cure conditions, etc.) under identical environments (temperature, 
moisture, etc.) and loading conditions (load rate, R-ratio, frequency, etc.) as 
the* structure. The influence of other material, environmental, and loading 
variables has been examined [20-22]. However, most of this work has been 
performed for static toughness and/or delamination growth. Much work still needs 
to be done to determine the influence of these variables on del ami nation onset 
in fatigue. 

Furthermore, the fatigue life methodology outlined in this paper for 
composite laminates took advantage of the real-time measurement of stiffness 
loss that was attainable in the laboratory for these simple coupon tests. For 
most structural applications, however, this information may not be available. 
However, in the absence of accurate damage growth and/or stiffness loss 
information, a similar, but more conservative, Damage Threshold/Fai L Safety 
approach may be used for composite structures [23]. 
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CONCLUSIONS 


A tension fatigue life prediction methodology for composite laminates is 
printed. Tension f.tlgu. tests were conducted on quas l- i sotropi c and 
orlhotropic glass epoxy, graphite epoxy, and gl ass/ graph i te epoxy hybrid 
aminates. Edge del am l nat l on onset data were used to generate plots of strain 
energy release rate as a function of cycles to delamination onset. These ploSs 
were then used along with strain energy release rate analyses of delaminations 
initiating at matrix cracks to predict local delamination onset. Stiffness loss 
was measured experimentally to account for the accumulation of matrix cracks and 
n ^iohT 310 " growth. Fatigue failure was predicted by comparing the increase 
in global strain resulting from stiffness loss to the decrease in laminate 

fami Ur \ fh" r ® 3ulting from delaminations forming at matrix cracks through the 

h ! 33 ' G °° d agreement between measured and predicted liv-s 

ndicated that the through-thickness damage accumulation model can accu-'tely 
describe fatigue failure for lamihates where the delaminatioo onset behavior in 

fatigue 13 well characterized, and stiffness loss can be monitored in real time 
to account for damage growth. 
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TABLE 1 ~ T190 E~Glass/X791 /50 Epoxy Lamina Properties 
(Batch //I) 

Avg. Properties Initial Reduced 


E n , GPa 

44.23 

44.23 

E 22 , GPa 

16.21 

11.50 

G 1? , GPa 

6.18 

3.09 

V 12 

0.284 

0.284 

e F 

0.0176 



TABLE 

2 - T 1 90 E-Glass/X751 /50 Epoxy Lamina Properties 


(Batch 

#2) 

Avg. Properties 

Initial 

E 1 1 ’ 

GPa 

39.44 

E 22’ 

GPa 

15.37 

G 12’ 

GPa 

5.10 

V 1 2 


0.270 

C F 


0.0256 
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TABLE 3 HM35 Graphite X7b1/50 Epoxy Lamina Properties 


Avg. Properties 

Ini tial 

Reduced 

E 1 1 , GPa 

152.50 

152.50 

E 22 , GPa 

7.50 

7.50 

G 12 , GPa 

4.^2 

2.06 

V 1 2 

0.246 

0.246 

e F 

0.00719 


TABLE 4 

- LAMINATE 

PROPERTIES 


MATERIAL 

LAYUP 

AVG. 

AVG. 

INITIAL MODULI 

MODULI BEFORE 

DELAMINATION 

e f 

E o 

(GPa) e lam 

E Q (GPa) 

e lam 

Glass 

[±45/0/90] 

s 

0.01500 

23.9 

25.2 

21 ,2 

21 .0 

Graphi te 

[±45/0/90] 

s 

0.00663 

57.5 

56.7 


— 

Hybrid 

[±45/0/90] 

s 

0.00695 

52.1 

51 .8 

— 

— 

Hybrid 

[±45/0] 

s 

0.02340 

24.4 

25.6 

21 .5 

21 .6 


20 


i 




21 


Tensile Stress-Strain Behavior of Graphite Epoxy and Glass Epoxy 
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Fig. 2. Tensile Stress-Strain Behavior of Graphite/Glass Epoxy Hybrids 




Fig. 3. Fatigue Damage in Glass Epoxy Quasi-isotropic Laminates 
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Fig. 4. Schematic of Delamination in a Quasi-isotropic Laminate Subjected to 
Cyclic Tension 
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Fig. 5. Tension Fatigue Behavior of (45/-45/0/90) s E-Glass Epoxy Laminates 
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Fig. 6. Tension Fatigue Behavior of (45/-45/0/90) s Graphite Epoxy Laminates 
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Tension Fatigue Behavior of (45/-45/0/90) s Graphite/Glass Epoxy 
Hybrid Laminates 
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Cycles, N 

Fig. 8. Tension Fatigue Behavior of (45/-45/0) s Graphite/Glass Epoxy Hybrid 
Laminates 
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Fig. 9. Micrograph showing Local Delamination From Matrix Crack in 
(45/-45/0) s Hybrid Laminate 
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Fig. 10. Influence of Damage on Laminate Stiffness 
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Fig. 11. Influence of Constituent Material on Stiffness Loss Associated with 
Matrix Cracking 
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Fig. 12. Stiffness loss due to delamination. 
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Fig. 13. Prediction of Stiffness Loss in Composite Laminates 
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90 degree ply Matrix Cracks 



= strain at failure in static test 
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Fig. 15. Increase in Global vs. Local Strain 
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ocal Strain Concentrations due to Through-the-thickness 
Accumulation of Local Delaminations in (45/-45/0) s Laminates 
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Fig. 17. K e Variation Across the Laminate Width in (45/-45/0/90) s Laminates 






Fig. 18. Local Strain Concentrations due to Through-the-thickness 

Accumulation of Local Delaminations in (45/-45/0/90) s Laminates 
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Fig. 19. Local Strain Concentrations in Zero Degree Plies due to Delamination in 
Quasi-isotropic Laminates 
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Fig. 21. Mechanical Strain Energy Release Rate at Delamination Onset as 
Function of Fatigue Cycles 
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train Energy Release Rate for Local Delamination Onset 
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Fig. 23. Reduction in Global ef due to Local Delamination Through the Thickness 
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Fig. 24. Tension Fai 
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Fig. 25. Delamination Onset Criterion for (45/-45/0/90)$ E-Glass X751/50 Epoxy 
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Fig. 26. Stiffness Loss as a Function of Delamination area for HM35 Graphite- 
X751/50 Epoxy 
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Fig. 27. Delamination Onset Criterion for (45/-45/0/90) s Graphite X751/50 Epoxy 
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Fig. 28. Comparison of Delamination Onset Criteria for E-Glass Epoxy and 
Graphite Epoxy 
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Fig. 29. Delamination Onset Criterion for (45/-45/0/90) s Graphite/E-Glass 
X751/50 Epoxy Hybrid Laminates 
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Fig. 30. Comparison of Delamination Onset Criteria for (45/-45/0/90) s E-Glass 
Epoxy and Hybrid Laminates 
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Fig. 31. Fatigue Life Prediction Based on Local Delamination Accumulation 
through the Thickness for (45/-45/0/90) s E-Glass Epoxy Laminates 
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Fig. 33. Fatigue Life Prediction for (45/-45/0) s Graphite/E-Glass Epoxy Hybrid 
Laminates 
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